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0. The aim of the present article is to reveal a structure shared by two basic zeta-functions in 
their fourth power moments. It might induce one to ponder over the possibility to go beyond. 

To be precise, let -T be a discrete subgroup of a Lie group G in the framework of GL2, 
and let be the L-function associated with a i^-automorphic function if) on G. Then a 
major subject in analytic number theory is offered by the mean value 

/oo 
\L^ {\ + it)\^ g{t)dt. (0.1) 
-00 

The principal issue is to establish an explicit spectral decomposition, and the following two 
cases have so far been discussed in greater detail: 

(1) the fourth power moment of the Riemann zeta-function ([9] [24] [25] [30]), 

(2) the fourth power moment of the Dedekind zeta-function of the Gaussian number field 

([8]), 

which correspond, respectively, to the specifications 

G' = PSL2(M), r = PSL2(Z), (0.2) 
G' = PSL2(C), r = PSL2(Z[v^]). (0.3) 

We shall see that in these cases the spectra come from irreducible representations of G oc- 
curring in LF'{r\G), and that the resulting integral transform of g has a kernel composed of 
Bessel functions of representations of G. It should be noted that (2) can readily be extended 
to any imaginary quadratic number field of class number one. Other examples that share the 
same structure and have been more or less worked out are 

(3) the mean square of the Dedekind zeta-function of any quadratic number field ([29]), 

(4) the fourth power moment of the Dedekind zeta-function of any real quadratic number 
field with class number one ([6]), 

(5) the spectral fourth power moment of all Hecke L-functions, under either (0.2) or (0.3) 

([34]), 

As far as our present purpose is concerned, (1) is the most fundamental, and (2) comes next 
endorsing our conceptual view about M(L^, (7) in general. The situation with (3) and (4) is 
much similar to (1), though (3) requires F be replaced by a Hecke congruence subgroup, and 
with (4) we need to move to Hilbert modular groups. The case (5) might appear different from 
others, but can in fact be regarded as an extension of either (1) or (2), though the situation 
with (0.3) is still under investigation. Note that in (2)-(4) as well as in (5) under (0.3) the 
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twist and the average with respect to Grossencharakteren can also be taken into account. 
Concerning this, an interesting argument has been developed by P. Sarnak [37], which seems 
to indicate another way to view (2)-(5). Despite this, we shall mostly concentrate on (1) and 
(2). 

To the list above one may wish to add, for instance, the mean square of individual 
Hecke L-functions associated with cusp forms, and further more any extension to the SL3 
environment ([10]). It appears, however, somewhat premature to discuss these subjects fully 
in the perspective described or suggested in the present work. We plan to return to them in 
the near future; nevertheless, see [18] [26] [31] [35] [36] for instance. 

Convention: The weight function g in (0.1) is assumed to be even, entire, real on M, and 
of rapid decay in any fixed horizontal strip. The symbol Q is reserved for the Riemann zeta- 
function, as usual. We shall often use F to denote various functions which are specified in 
local contexts. Other notations are introduced where they are needed for the first time, and 
will remain effective thereafter unless otherwise stated. Also we stress that we are concerned 
with the structural aspect, and the asymptotic study is being more or less left aside; by the 
same token we shall often skip the discussion of convergence, naturally under the premise 
that no confusion be brought in. 

Acknowledgement: This article is an outcome of our recent works which we conducted 
either solely or jointly with R.W. Bruggeman. We are greatly indebted to him for his invalu- 
able cooperation. The text below will thus contain various excerpts from those works, with 
appropriate modifications. We thank A. Ivic and M. Jutila for their constant encouragement. 

1. First we shall describe an old observation ((2.3)-(2.5) below) made in [24] (see also [30, 
Section 4.2]), a proper exploitation of which was only recently actuated in our joint work [9] 
with Bruggeman on M(C^,^) the fourth power moment of the Riemann zeta- function. 

Thus we begin with an idea of H. Weyl: In analytic number theory we deal mainly with 
sums over rational integers 

E^H' (1-1) 

where / is an interval in M, and F an arbitrary function; and to estimate this there is a 
general principle due to Weyl. A version of it is attributed to J.G. van der Corput and built 
upon the following triviality: 

^ 00 M 
ne/ n= — oo m=l 

where M > 1 is arbitrary, and 5j the characteristic function of /. If some effective inequalities 
are applied to the right side, then the result does not look trivial any more but can be a sharp 
tool (see [16, Chapter 2]). In fact it could yield subconvexity bounds for Q which have played 
basic roles in many problems in analytic number theory. For instance the bound 

C(| + it) <t^/^logt, t>2, (1.3) 

is relatively easy to prove with (1.2), which is significantly better than the convexity bound 
following from the functional equation (1.7) below. 



Mean Values of Zeta-Functions via Representation Theory 3 

One may regard (1.2) as a kind of lifting of a one dimensional sum to a two dimensional 
sum. The original problem has been transformed into the one of finding interaction among 
the non-diagonal entries. This observation leads us to another triviality: For general double 
sums we have the decomposition 

^ F{m, = { E + E + E }^(^' (1-4) 

m,n m=n m<n m>n 

This is exactly the same as what F.V. Atkinson did in his important investigation [1] on the 
mean square of In [30, Section 4.1], an essentially equivalent assertion is formulated as 

M(C,^) = j Re |^(i+it)|+27£;-log(27r) g{t)dt + 2T,Re [g 

+ 4y'c/(n) / (r(r + l))"^/V(log(l + l/r))cos(27rnr)dr, (1.5) 

n=l 

where 7_b is the Euler constant, d{n) the number of divisors of n, and 

/oo 
g{t) cos{xt)dt. (1.6) 
-oo 

We should note that (1.5) implies (1.3); see Ivic's lecture notes [15] for the details of various 
consequences of Atkinson's result. We stress also that the last infinite sum, which is of 
a spectral nature, corresponds to the non-diagonal parts in (1.4), and that (1.5) could be 
regraded as a completion of a particular application of Weyl's idea to 

This completion was realized in [30] via the functional equation 

C(l -s) = 2(27r)-^ cos (Ins) r(s)C(s). (1.7) 
On the other hand, in [1] had been applied the Poisson summation formula 

oo «(X) °° roo 

VF(n)=/ F(r)dr + 2V / F(r) cos(27rnr)dr, (1.8) 
n=i -^0 n=i -^0 

where F is assumed to be smooth and of fast decay on (0, cxd). However, this difference 
is superficial, for (1.7) and (1.8) are equivalent to each other. In fact, the kernel function 
cos(27rr) in (1.8), as well as in (1.5), is related to (1.7) via the Mellin transform: 

/■oo 

/ cos(27rr)r*-Mr = (27r)-*cos (|7rs)r(s), 0<Res<l, (1.9) 
Jo 

or, more precisely, via the convolution relation 

/ F(r) cos(27rnr)dr = ^ I F*{l-s) cos (^tts) T{s){2TTn)~^ds, (1.10) 

where F* is the Mellin transform of F, and {a) the vertical line Res = a with o; > 0. 
Summing (1.10) over n, and applying (1.7), we are led to (1.8). Reversing the reasoning, one 
may reach (1.7) with a combination of (1.8) and (1.9). 
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In terms of the representation theory of the Lie group R, the identity (1.8) is understood 
as a transformation formula of F that gives a way to compute the values of projections of the 
Poincare series 

YF{n + x), xeM, (1.11) 

to irreducible subspaces of L^(Z\M). The appearance in (1.8) of the Fourier-cosine transfor- 
mation indicates the nature of the harmonic analysis over the group M in which Z is a discrete 
subgroup. 

A motivation of the present article is to see how far extends the above harmonic structure 
lying behind the important expansion (1.5). 

2. In this context, one may ponder whether the dissection mode in (1.4) is optimal or not. It 
is certainly not in general. An option to tune it up is to refine the notion of being diagonal. 
A simple arithmetic way to do this is to replace (1.4) by 

J]F(m,n) = | E + E + E \nm,n), (2.1) 

m,n \km=ln km<ln km>ln) 

where A;, / are arbitrary non-zero integers. To extract information from all of these dissections 
we multiply both sides by a weight W{k,l) and sum over all k, I. We get 

J2w{k,l) \ (E^K"))=| E + E + E \w{k,l)F{m,n). (2.2) 

. k,l / \m,n J \km=ln km<ln km>ln) 

As before, we may regard (2.2) to be a lifting of a double sum to a four dimensional sum; 
and we are led to the following trivial decomposition: 

E F{k,l,m,n)^l E + E + E \F{k,l,m,n). (2.3) 

k,l,m,n \km=ln kmKln km>ln) 

Then we take a new viewpoint: We regard quadruple sums as sums over 2x2 integral 
matrices M. The last identity thus becomes 

Y,nM) = \ E + E + E \nM). (2.4) 

M [|M|=0 |M|>0 |M|<oJ 

Invoking Hecke's representatives of integral matrices with a given determinant, we have further 

oo 

E nM) = EV7^/^i^/(l)' PF(g) = E^(7g)> geC, (2.5) 
|M|>0 /=! i^r 



with (0.2), where F/(g) = F{^/f ■ g) and 

(2.6) 



T,F(g) = -l=EE^ 
d\f b=i 



y7/d b 
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Hereby we find a relation lying between the dissection (1.4) and the theory of r'-auto- 
morphic functions on G. If F is sufficiently smooth one may apply the harmonic analysis over 
r\G to the Poincare series Pp- Thus, we are now to find the transformation formula of F that 
gives a way to compute the values of projections of Pp to irreducible subspaces of L'^{r\G). 
In other words, our interest is in fixing the kernel function corresponding to cos(27rr) of (1.8), 
via the theory of irreducible representations of G occurring in L'^{r\G). In the application 
to M(C^,5'), as developed in [9], those four integral variables in (2.3) correspond to the four 
zeta- values in an obvious way. In the subsequent sections we shall show the salient points of 
the pertinent parts of [9] and briefiy describe a solution to the last problem. 

It should be stressed here that the restriction of the sum (2.5) to those M with \M\ = /, 
i.e., TfPpfi is perhaps more worth investigating. A typical example is the additive divisor 
sum 



Y,CTx{n)a^{n + f)W{n/f), 



(2.7) 



where (7\{n) = (i^ with A e C; and the smooth weight W has a support in the positive 
reals. A detailed discussion of the spectral decomposition of (2.7) is developed in [27], and 
the results there have recently found important applications in our joint work [19] (see also 
[20]) with M. Jutila, where is considered an instance of the case (5) above, and established 
a new uniform bound for Hecke L-functions associated with cusp forms under (0.2), which 
generalizes the classical subconvexity bound (1.3) to a vast family of L-functions. 

3. In this and the next two sections, we shall deal with the Poincare series Pp, with an analogy 
of (1.8) in mind. To begin with, we shall collect elements of the theory of T-automorphic 
representations of G, under the specification (0.2): Thus we write 



n a; 



1 



a[y] 



k[^] 



cos 9 
— sva.9 



sin 9 
cos 9 



(3.1) 



Let N = {n[x] :xeR},A = {a[y] :y>0}, and K = {k[9] : 9 e K/ttZ}, so that G = NAK 
be the Iwasawa decomposition of G. We read it as G 9 g = nak = n[a;]a[y]k[^]. The Haar 
measures on the groups A^, A, AT, G are defined, respectively, by dn = dx, da = dy/y, dk = 
d9/7r, dg — dndadk/y, with Lebesgue measures dx, dy, d9. The space L'^{r\G) is composed 
of all left T-automorphic functions on G, vectors for short, which are square integrable over 
r\G against dg. Elements of G act unitarily on vectors from the right, and we have the 
orthogonal decomposition into invariant subspaces 

(r\G) = c • 1 '^L^ (r\G) %^ {r\G) , (3.2) 

where is the cuspidal subspace, and '^L'^ is spanned by integrals of Eisenstein series. 
The cuspidal subspace splits into irreducible subspaces: 



(3.3) 



where = (5^ + dy) — yd^de is the Casimir operator. Under (0.2), we can restrict our 
attention to two cases: either vy e z(0, oo) or vy is equal to half a positive odd integer. 
According to the right action of the space V is decomposed into A-irreducible subspaces 



y= 



dimFp < 1. 



(3.4) 
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If it is not trivial, Vp is spanned by a i^-automorphic function (pp such that (/:p(gk[6']) = 
exp{2ip9)(p{g); it is called a r'-automorphic form of spectral parameter uv and weight 2p. 

Let us assume temporarily that V belongs to the unitary principal series, i.e., G 
i{0, oo) under (0.2). Then dimT^ = 1 for all p e Z, and there exists a complete orthonormal 
system {(pp eVp-.peZjoiV such that 

oo 

<^p(g) = ^ \n\~''gv{n)An(l)p{g;i^v) 



n= — oo 
n/0 



n= — oo 



= Yl ^^/T-T •^sgn(n)</'p(a[|n|]g; uv), (3.5) 

'\n' 



where 0p(g; i') = y^/'^'^'^ ex.p{2ip9), and Au is the Jacquet operator: 

Au4>p{g,i') = / exp{—27riuv)(j)p{wn.[v]g;i')dii, w = k [^tt] . (3-6) 
Jn 

It should be observed that the coefficients Qvin) in (3.5) do not depend on the weight. We 
note that for u e 

AuMs-^u) = y'/'-^eM2^iux)eM2piO) ^^^^ (^)'^^ 

= (_l)p^i/2+^|^|'-i/2exp(27rmx)exp(2pi^) ^^Mu)pAMu\y) ^3 7) 

r(sgn('u)p+ ^ + 1^) 

where W\^^{y) is the Whittaker function (see [WW, Chapter XVI]). The first line is valid for 
Hjev > 0, while the second defines An(f>p for all v & C In particular, we have the expansion 

¥'o(g) = 7^7^ ^^/y V Qv{n)K^^{27r\n\y)ex.p{27rinx), (3.8) 

^ ^ n— — oo 

with being the i^-Bessel function of order u. This is a cusp-form on the hyperbolic upper 
half plane G/K. 

Next, let us consider a F in the discrete series; that is, — ^ — ^, 1 < £ G Z. We have, 
in place of (3.4), 

~oo -e 
either V = ^Vp or V= ^ Vp, (3.9) 

p=£ p= — oo 

with dimV^ = 1, corresponding to the holomorphic and the antiholomorphic discrete series. 
The involution g = nak 1-^ n~^ak~^ maps one to the other. In the holomorphic case, we have 
a complete orthonormal system {(pp : p >£} in V such that 

np-e + i) ) Il^^~''^Mn)AnMs;^v)- (3-10) 
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In particular, we have 

92^ 1/2+^ 00 

^i{g) = (-1)^— ==exp(2z£^)/ J]^,.(n)n^-V2exp(27rm(x + zy)). (3.11) 

V ^ (^^) n=l 

This infinite sum is a holomorphic cusp-form of weight 2£ on G/K. 

With this, we may assume further that all V be Hecke invariant, so that there exists, for 
any integer n > 1, a real number Tv{n) such that 

Tn\v = Mn) ■ 1, (3.12) 

where T„ is as in (2.6). Thus, for any non-zero integer n, 

Qv{n) = Qv{sgn{n))Tv{\n\). (3.13) 

We may introduce the convention that qv{—^) = and ^v'(l) = for F in the holomorphic 
and antiholomorphic discrete series, respectively, and gv{—^) = ^vQv{^) with e^/ = ±1 for 
V in the unitary principal series. We associate, with each V, the Hecke series 

00 

Hv{s) = J2^v{n)n-', Res > 1, (3.14) 



n=l 



which continues to an entire function. 

4. With this, we are going to decompose Pp via (3.2). We may restrict ourselves to the 
orthogonal projection wy with V in the unitary principal series. In fact, the discrete series 
is highly analogous, the projection to the subspace ^L^(r\G) is facilitated by that Eisenstein 
series are explicitly defined, and the space of constant functions gives no specific problem. 

We are of course concerned with how to compute zjvPf in terms of F, or more generally 
with the harmonic analysis on V while the variable g e G is restricted to the big Bruhat cell. 
To this end we shall employ a reasoning which we term the Kirillov scheme. This is because 
we utilize an operator %, defined by (4.3) below, that apparently originated in A. A. Kirillov 
[21]. We shall exploit two basic properties of X, and they are embodied here in two lemmas, 
respectively: 

Lemma 1. Let U = U^, with 1/ & iM. be the Hilbert space 



p= — oo 

equipped with the norm 



u 



|cpp, 5^ Cpcl^p- (4.2) 

\p=- 



p=— 00 p=— 00 
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For u e and smooth (j) E U, that is, with Cp decaying faster than any negative power of 
\p\, we let 

X(l>{u) = \u\^/^-''A^(l){l) = yisgn(n)</'(a[|w|]). (4.3) 
Then the operator % maps U unitarily onto L^(R^, {l/7r)d^), where d^u = du/\u\. 
Lemma 2. Let us define the Bessel function of representations o/PSL2(lR) as to he 

Uu) = (J-Ti"^(47rv1^) - J,T^"^ (47^^M)) , (4-4) 

with Jjf = Ji, and J~ = in the ordinary notation for Bessel functions. Then, for any 
smooth (f) & Ui,, we have 

%R^4>{u)= I j^{uv)X(f){v)d''v, weM^. (4.5) 



Proof. It should be stressed that the definition (4.3) is taken from [9] [35], and somewhat 
different from that employed in [32] [33]. A proof of the unitarity of % is given in [32] [33, 
Theorem 1]. It depends on the following integral formula via the second line of (3.7): For any 
a,PeC and |Rez/| < | 



TT 



Jq u [X — jji) sm(27rz^) 

1 

X 



1 



(4.6) 



which is tabulated as [14, eq. 7.611(3)]. The verification of this made in [32] [33] employs 
the Whittaker differential equation ([41, p. 337]) which is related to the Casimir operator. 
The surjectivity of % is proved in [9] via the first line of (3.7) and the completeness of the 
system {{{v + i)/{v — i)Y : p G Z} in the space (M, dv/{7r{v'^ + 1))). As to Lemma 2, the 
realization (4.5) of the action of w the Weyl element in terms of the space L^(R^, {l/7r)d^) 
seems to have been published for the first time by N.Ja. Vilenkin (see [38, Section 7 of 
Chapter VII] as well as [39, eq. (17) on p. 454]), though the concept of the Bessel function of 
representations had been coined by I.M. Gel'fand, M.I. Graev and I.I. Pyatetski-Shapiro [12]. 
Two independent proofs are known; they are conceptually different. One is due to M. Baruch 
and Z. Mao [2], which is along the line of [38] and fills a gap therein concerning a convergence 
issue. The other is due to ourselves [31] [32] [33], and seems more in line with the purpose of 
the present work. It is shown there that (4.5) is in fact equivalent to the Jacquet-Langlands 
local functional equation ([17, Theorem 5.15]) 

{-l)PTp{s) =2i-2«7r-2*r(s + i^)T{s - i^) 

X (cos(7rs)rp(l -s) + cos(7rzv)r_p(l - s)) (4.7) 

for the Mellin transform 

/•oo 

rp(s)= / AuM^)u"~"~^du, (4.8) 
Jo 
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which continues meromorphically to C. The Melhn inversion of (4.7) coupled with (4.9)- 
(4.10) below gives (4.5) for (j) = (j)p; the extension to any smooth e C/ is easy. Thus, if 
|Rez^| — I < Res, then 



/. 



° 1 



j^{u)\u\'-^du = -(27r)-2^ cos(7riy)r {s + ^ + u) F {s + ^ - u) ■ (4.9) 



and if |Rej^| — | < Res < —j, then 



r 1 

/ j^{u)u'-^du= (27r)-2"sin(7rs)r(s+ i + i^)r (s + i -zy) ; (4.10) 

Jo 

The former follows from [40, eq. (8) in Section 13.21], and the latter from [ibid, eq. (1) in 
Section 13.24]. This ends our brief discussion on the proof of the above lemmas. 

The above extends not only to the discrete and the complementary series but also to the 
complex situation, i.e., to PSL2(C), as is to be shown in Section 7. 

5. Now we shall carry out the computation of wvPf via the Kirillov scheme. It should be 
stressed that absolute convergence required below can readily be confirmed, provided F is 
sufiiciently smooth. 

Thus, the projection to Vp is, by the unfolding argument, 

{PF,'fip)r\G^ / F{g)ifp{g)dg 

JG 

where (3.5), (3.13) are used; Fm(g) — F{a,[m\~^g) and 

^lF{v) = f F{g)AsMs)dg. (5.2) 
Jg 



Thus 



p\S) 



p=— oo 

oo oo 

= Mir EE 



Tv{m)Tv{n) 



where 



X (^^(+^+) + + ev^S(+'-) + evS(-'+)) F^(a[n]g; lyy), (5.3) 

oo 
p= — oo 

oo 

= eM'^TTi62x) J2 ^pF{u)As,(l)p{a[y])exp{2ipe). (5.4) 



p=— oo 
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Since it can be asserted, with an appropriate change of F, that our interest is in the 
value wvPf{^) (see (2.5)), we may restrict ourselves to the subgroup A. Namely, it suffices 
to consider 'B'^^'^'^^^ F{a[y\] v), with a new F; and this can be expressed in terms of the Kirillov 
operator: 

oo 

'h^'-^'^^F{4y]-v) = %L'-F{52y), F = J] $jF(z.)0p, (5.5) 

p= — oo 

where L^F is smooth. Hence, we may proceed in the sense of weak convergence, appealing 
to Lemma 1; and we observe 

^lF{v) = {Ju^F,(f>p)u =- [ XL^F{u)X^p{u)d''u. (5.6) 

This means that if we are able to transform (5.2) into 

^lF{v) = - [ Y\u)X(l)p{u)d''u, (5.7) 

then it should follow that 

'B('^''^^F{#,iy) = Y'^{S2y), (5.8) 

because of the surjectivity assertion in the lemma. 

Since the integral in (5.2) is in fact over the big Bruhat cell, we perform the change of 
variables accordingly. We have instead 

^'pF{u) = / / F{a[u]g)R^AsM^H)dg—- (5.9) 
Jo JNwN 

Here i?g is the right translation with g = n[a;i]wn[a;2], and dg — dxidx2/iT. We observe that 

RgAs(f)p{a[u\) = ex.p{2'Ki5xiu)AsRwRn[x2]4'p{^b^]) 

= exp{27riSxiu)%RyfRn[x2](l^pi^'^)- (5.10) 

By Lemma 2 this is replaced by 

RgAs(f)p{a'[u\) = exp{27ri5xiu) / exp{27rix2v)j,^{5uv)%(f)p{v)d\, (5-11) 

Jrx 

and (5.9) by 

^^F{u) —- / / F{a.[u]ii[xi]wn[x2])exp{-2m5xiu) 

X / e'x.p{—2-nix2v)jy{5uv)%(j)p{v)d^vdxidx2 — • (5-12) 

Hence we find via (5.8) that 

^(5i,5.)^(a[y];iv)= / U5^52yu) 
Jo 

x \ [ F{aL[u]n[xi]wn[x2]) exp{—2ni6iuxi — 2niS2yx2)dxidx2\ —, (5.13) 

./»2 u 
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which ends the apphcation of the Kirillov scheme. 

We may compare (4.7) with (1.7). Then (4.5) may also be compared with (1.8). That is, 
the formula (5.3) coupled with (5.13) which is a local assertion derived from (4.5) corresponds 
to (1.8). As remarked above already, analogues of Lemmas 1 and 2 are shown in [9] [32] for the 
discrete series representations and the complementary series, although the latter is irrelevant 
under (0.2). Hence, returning to (3.2), we obtain a genuine extension of (1.8). The final result 
is, however, too complicated to be stated as an independent assertion. We should instead 
be content with the local expression (5.13) and with the fact that we have found that the 
combination of Lemmas 1 and 2 is the key implement. 

Hence, what corresponds to the cosine-transform in (1.8) is (5.3) with (5.13). One may 
desire to compute the double sum (5.3) and the last double integral into closed forms. In 
the applications to M(C^,5') and to the sum (2.7), which will be briefiy dealt with below, 
we are in a fortuitous situation that the double sum is transformed into a product of two 
values of Hy- As to the double integral, it is a Fourier transform over the Euclidean plane, 
and thus might be expressed in terms of a Bessel transform. With M(C^, (7) as well as (2.7), 
the situation turns out in fact to be as such, and we shall see that (5.13) is expressed as an 
integral transform whose kernel is a convolution of two instances of the Bessel function of 
representations (see (6.2) below). 

Thus the matter seems to depend much on the specific nature of the seed F. Nevertheless, 
with any smooth F, one might appeal to Mellin transform of several variables, and the above 
could be pushed into a more closed form. 

6. We are now at the stage to render the spectral decomposition of M(C^,(7) in terms of 
notions from representation theory: Thus, let us put 

^^-^^^ = J^)L (^) ^e(log(l + l/.))=(.;.)dx., (6.1) 

S(«;^)= / jo{-v)jJ-)^. (6.2) 
Then we have 

M{C\g) = {m(^) +M(=) +M(^)} (C^^), (6.3) 

where 

M^'^HC',^) = {^fe{uv;g), (6.4) 

V 

with ay = |^\/(1)P + |^'i/(— 1)1^- The V runs over a maximal orthogonal system of Hecke- 
invariant cuspidal F-automorphic representations of G. Apart the term 27rRc {(log(27r) — 
7£;)5'(|^) — l^fi^'dO}' '^^^\C'^i9) is an integral transform of g whose kernel is given 
explicitly in terms of logarithmic derivatives of the Gamma function. 

The proof of (6.3) as developed in [9] starts with the integration of C{zi + it)C{z2 + 
it)({z3 — it)({z4^ — it) against g{t)dt over R, where (2^1, 2:2, -^3, -^4) is to remain in the region 
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of absolute convergence. The device (2.3)-(2.5) is to be applied. However, the naive choice 
of the seed does not work well, because we require it be smooth and of rapid decay on G. 
Such a choice of the seed is not difficult but somewhat subtle. Thus, we are forced to employ 
instead a sequence of suitable F's and a limiting procedure with respect to F. Nevertheless, 
those F chosen in [9] are nice in that they are ^-equi variant, i.e., F(a[y]g) = y'^F(g) with 
an a; G C being independent of F. Hence the summation over / in (2.5) is the same as the 
multiplication by a value of Hy at each V; and also the sum over m in (5.1) can be written in 
terms of a product of a value of Hy and ($+ + ey$~)F(i^y). In particular, the omission of m 
in (5.4) is possible without changing the specification of F, which amounts to a considerable 
simplification in the subsequent discussion leading to (5.13). Besides, this makes (5.13) easier 
to handle; taking the limit in F we come already close to the expression (6.2). Moreover, the 
sum over n in (5.3) yields now another factor in a value of Hy- Then it remains to perform 
analytic continuation and specialization with respect to {zi, Z2,zs, Z4). Thus the factor Jq in 
(6.2) stands for the integral kernel of the Bessel transform that emerges from the inner double 
integral of (5.13) at the end of the whole procedure (see [9, (7.34)]). This explains how (6.2) 
originates and reveals especially the mechanism behind (6.4). 

The same argument can be applied to the additive divisor sum (2.7). We are led to a 
spectral decomposition analogous to (6.3). Specifically, the cuspidal part is found to be 

l/(A+p+i)/2 J2 avry{f)Hy (i(l - A - /.)) i^v (1(1 + A - /.)) 

V 

x{^++ey^-){uy;X,ii;W), (6.6) 

where 

/■oo 

^s{i';\l^;W)= / W{u)As{u;i^;X,iJ,)u^^+'^^/^+^d''u, (6.7) 
Jo 

with 

A5{u;iy;X,iJ,) = j\/2{-5v)j^,{5v/u)-^—^. (6.8) 

It is safe to keep both |ReA|, |Re//| sufficiently small so that (6.6) holds with the expression 
(6.8). However, as can be seen from (4.9)-(4.10), can be expressed in terms of the Mellin 
inversion of a product of four Gamma factors, and then (6.7) allows us to continue (6.6) 
analytically to quite a wide domain of (A,//). 

It has been explained in the above how the factor jq in (6.2) turns up; the same can be 
applied to j\/2 in (6.8). However, it is not done in any framework of metric theory. This is 
sharply different from the situation with another factor jj^ shared in both equations. Thus 
it remains still to find a genuine characterization of the factors jo, j\/2- In passing, we note 
that instead of (6.2) we may write 

S(u; .) = 2R. ju-/-" (1 - -J-^) Jfe^^""' + + 2.; -l/„) |. (6.9) 

with the Gaussian hypergeometric function 2F1 (see [30, (4.7.2)]). This reminds us of the free- 
space resolvent kernel of the hyperbolic Laplacian (see [ibid, (1.1.49)]), a fact that appears 
mysterious to us. 
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It might be expedient to make here a digression on a historical background: A prototype 
of the spectral decomposition of M((^^, g) was obtained by the present author in [24] [25], which 
was afterwards improved to (6.3) in [30, Theorem 4.1]. However, the assertion there did not 
reach the expression (6.2); it was stated with (6.9). A reason for this is in that there we used 
the Kloosterman-Spectral sum formula of N.V. Kuznetsov [30, Theorems 2.3 and 2.5], which 
is pretty handy but hides the mechanism working behind the integral transform appearing on 
the spectral side. Note that in the above we dispensed with Kuznetsov's sum formula. The 
argument of [9], whose most salient part is depicted in the previous section, is admittedly 
more involved than that in [30] , but this is much due to the fact that we started from the very 
fundamental assertion (3.2), whereas the discussion in [30] lacks the perspective offered by 
representation theory. Thus, the appearance in Kuznetsov's sum formula and consequently 
in [30, Theorem 4.1] of the contribution of holomorphic cusp forms was just an accidental 
byproduct of a technical marvel and remained mysterious there. Our discussion of M{('^,g) 
in terms of a Poincare series on the group G allows us to see all contributions of cusp forms in 
a fairly equal term, since our method is based on (3.2), where all irreducible representations 
have equal rights. 

It was Bniggcman [3] [4] who tried for the first time to understand, via (3.2), all the 
terms on the spectral side in Kuznetsov's sum formula. However, the real comprehension of 
the structure supporting the sum formula appears to have been done by J.W. Cogdell and 

1. Piatetski-Shapiro in [11]. In particular, the Kirillov scheme together with the role of the 
Bessel function of representations was developed there, and Kuznetosv's sum formula was 
newly proved, though their discussion appears sketchy to us. The authors of [9] were inspired 
by the work [11]. 

7. The aim of this and the next sections is to show that the above discussion extends to the 
situation (0.3). In particular, we are going to show the complex analogues of Lemmas 1 and 

2. Note that some symbols used under (0.2) are now assigned to corresponding notions under 
(0.3); this convention should not cause any confusion. 

Thus, let G = PSL2(C), and put 



n[z] 



1 z 
1 



, h[..] 



u 



1/u 



a (3 
—P a 



(7.1) 



where z,u,a,P G C with u ^ 0, |ap + = 1; and also N — {n.[z] : z e C}, A — {a[r] : 
r > 0}, K = PSU(2) = {k[a,P] : e C} with a[r] = h[y/r]. In terms of the Euler angles 
(fi, 9, i/j, we have 



k[a, P] = h[e'^/^]v[i9]h[e''^/% y[e] = k[cosh(^/2), sinh(^/2)]. 



(7.2) 



The Iwasawa decomposition G = NAK is read as G 9 g = n[2;]a[r]k[a, (3]. The Haar measures 
on respective groups are given by (in = dz^ da, — dr/r, dk = sin 9 d(pd9 dip /{8tt^), and dg = 
d-n.ds.dk /r"^. With this, the Hilbert space L^{r\G) is formed, to which G acts from the right; 
and we have an exact analogue of (3.2). See [5] [8] for more details of what follows. 

The cuspidal subspace decomposes into irreducible subspaces 



%\r\G) = ^v. 



(7.3) 
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To classify representations V, we need two Casimir operators Q±, = Q+, where 

11 1 rp'^'^ 

2 2 2 2sin^ 

+ \r^dl - \xrdrd^ - ^-dl - ^-rdr + \id^. (7.4) 
They become constant multiphcations in each V: 



^^±|v=xj •l = ^((pyTz^i/)'-l)-l, e Z, e «[0,oo). (7.5) 

The pair (py , uy) is called the spectral parameter of F, but in the sequel we shall write simply 
{pv.i^v) = {p, i^)- 

According to the action of the space V decomposes into K-irreducible subspaces 



y= yi,,^ dimy^,, = l. (7.6) 

\p\<l,\q\<l 

To describe this precisely, let VLk be the Casimir element of the universal enveloping algebra 
of K defined by 

— ^ idl + sin^ Odl + 91-2 cos 9d^d^ + sin 9 cos 9 de) . (7.7) 
2sin ^ ^ 

Then 

Vi,q = {FeV: QkF = + 1), d^F = -iqF} . (7.8) 

Any non-zero element of Vi^g is called a r'-automorphic form of spectral parameter (p, u) and 
K-type {l,q). 

Next, we define functions q on K by 

The system {^p^q ■ \p\, \q\ < 1 < is a complete orthogonal basis of L'^{K) with norms 

1 / 2/ \^/V 2/ ^ "^/^ 



2 



»t,jk=^^,:j . (7.10) 



Here are some of its properties which we shall need later: Under the convention that $p ^ = 
if the condition |g| < I is violated, we have 

$^,^(k[a,/3]) = e-^^^-^^^^lq{^,[i9]), (7.11) 



= i{l - ?X,+i(v[z^]) + i{l + ?X,_i(v[z^]), (7.12) 
2^^|f^*L(^[^^]) = ^(^ - - ^(^ + (7.13) 

2^^-^<.(v[^^]) = ^(^ + P+ - ^(^ - P+ lXi,,(v[z^]). (7.14) 



2de^lq{^[i9]) = ^(/ + P + l)$J,+i,^(v[i^]) + i(/ - p+ l)$;_i,,(v[i^]) 
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For a verification of (7.12)-(7.14) see [5, Lemma 5]. 

With this, we put 0/ ^(g; v) = r^'^^^p q(k)/||$p ^Wk, and its Jacquet transform is defined 

by 

'^u(l>l,qis;j^) = exp{-27riB^{uv))(f)i^g{wn[v]g;iy)dn. (7.15) 
Jn 

Let (fii^q be a generating vector of Vi^q, so that {(pi,q '■ \p\ < I, \q\ < /} is a complete orthonormal 
system in V. Then we have 

ViAs) = ^ \i^\~''Qvioj)^uj(f>l,qiSi i^), ^ e ^W, (7.16) 

which is precisely an analogue of the first line of (3.5). 

Further, we define the Hecke operator for each non-zero f E by 

I"' I d\f 6 mod d / V J J / 

and we assume that all V are Hecke invariant so that there exists a real number Tv{f) such 
that Tf\V — Tvif) • 1. In particular we have, for all non-zero n e Z[i], 

gy{n) = gvil)in/\n\rTv{n), (7.18) 

We have 

rv{—n) = Tvin), Tv{in) = eyrvin), ty = ±1 (7-19) 
The Hecke L-function of the space V is defined by 

iiv{s) ^\Y^Ty{n)\n\-^\ (7.20) 

which continues to an entire function; note that Hy = whenever ey = —1. 

Now, the complex analogues of Lemmas 1 and 2 are as follows: 
Lemma 3. Let U = Up^i, be the Hilbert space 



Ccf>l,q, Cf>l,q{g) = cf>l,q{g;u) (7.21) 
\p\<l,\q\<l 



equipped with the norm 



u — 



\p\<l,\q\<l \p\<l,\q\<l 



For u E and smooth (f) E U, we let 

3C0(w) = \u\^-''{u/\u\)PAu(l){l). (7.23) 
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Then the operator % maps U unitarily onto L^(C^, {2/7r)d^), where d^u = du/\u\'^. 
Lemma 4. Let us define the Bessel function of representations o/PSL2(C) as to he 

= 27r^^-- {Jp^u{'^7^u)J-p-u{'27ru) - J-p+j,{2'Ku)Jp+u{'2'Ku)) . (7.24) 

Then we have, for any smooth cj) e Up^v, 

3Ci?w0(w^)= / jp,r,{uv)X^{v'^)d''v, weC^. (7.25) 

Here Jp-i,{u)J-p-u{u) is understood to be equal to (it/|tt|)^|it|~^^ J*_j^(m) J*p_j^(tt) where 
J*(tt) is the entire function that coincides with u~^Jfj,{u) when tt > 0. 

With this, the analogues of (6.1)-(6.5) for the Dedekind zeta-function Cfe of the Gaussian 
number field k = Q{i) can be rendered as follows: Let us put 

Q(p,iy;g) = -—^^^ I y^^g^{2log\l + l/u\)E{u;p,u)d''u, (7.26) 
16sin(7rz/) l^t + 1| 

E{u;p,u) = f jo,o {y/^)jp,i^ (Vv/u) (7.27) 

Jcx \ / \V\ 

Then we have 

mlg) = +M(^) +M(^)} iClg), (7.28) 

with V running over a maximal orthogonal system of Hecke-invariant cuspidal P-automorphic 
representations of G. Here 'M^'^\Ck^ 9) is analogous to M^'^\C'^,g), and 

^^'\Cl 9) = Y1 \QvW\'Hv i^f 0{PV, i^V, 9), (7.29) 

V 

p^ooh^) |a(l + z^,2p)| 27rz 

where Ck{s-,p) is defined by 

Ck{s,p) = \ E(n/|n|)^P|n|-2^ Res > 1, (7.31) 

which continues meromorphically to C. 

The spectral decomposition (7.28) was proved in [8]. The argument was a faithful exten- 
sion of the older proof of (6.3); that is, it depended on the sum formula of Kloosterman sums 
under the situation (0.3) that was established also in [8]. Thus, in much the same mechanism 
as Lemmas 1 and 2 did with Kuznetsov's sum formula, the last two lemmas should allow us 
to dispense with the sum formula of Kloosterman sums, in deriving (7.28). This claim is still 
to be checked fully; but it is certain such a new proof is available and conceptually simpler 
than that in [8] . 
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Comparing (7.26)-(7.30) with (6.1)-(6.5), the outward similarity is striking. However, 
our good luck ends there. That is, the asymptotic nature of (6.3) is much superior than that 
of its counterpart (7.28). In fact, (7.28) does not seem suitable to be utilized as a means to 
derive quantitative assertions on the fourth power moment of Ck- Concerning this difficulty, 
Sarnak kindly suggested us to try to take further averaging: 

oo 

J2 nCl{;Q),9)h{q), (7.32) 

q=—oo 

with a smooth weight h. The argument of [8] should extend to this sum. The same can be 
said about an obvious analogue of (2.7). The latter is expected to yield an extension of the 
main result of [19] to L-functions associated with the group PSL2(Z[z]). In fact this has been 
in our current investigation. However, we have not worked out all the details yet. 

8. Now, we are about to prove Lemmas 3 and 4. This section is a reworking of our joint work 
[7] with Bruggeman; a few corrections are made. It should be stressed that the surjectivity 
assertion in Lemma 3 is a new addition, and that the definition (7.23) is somewhat different 
from that employed in [32]. 

We begin with the unitarity of %: Naturally it is sufficient to show the orthogonality 
relation: 

-/ X4>i,MX^ir^d''u = Si,rSq,q^, (8.1) 
with Kronecker deltas. By definition, 

Xcf>iju) = iu/\u\)-'^A,cf>i,MM)/\\^ljK. (8.2) 

Note 

•/li0«,g(g) = exp(27rzRe(2;)) ^m(0*L,g(k) 

\m\<l 

= exp(27rzRe (z)) e-^^'^-^^^^L(0<,,(v[^^]), (8.3) 

\m\<l 



where (7.11) has been used, and 
=-^i0/,m(a[r]) 



f exp(-27rzrRei;) , 
ic (1 + 1^1^)1+'^ 



-1 



V -1 



d+v 



(l + 7;2)i+- 



p,m 



vdv, (8.4) 



where d'^v = {dRev){dImv). The left side of (8.1) is equal to 

\\^p,q\\K Jo " 
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Then we observe that the functions Vq(r) satisfy the differential equations 

D-v[{r) = 47ri(/ + q)r-\[_^{r), (8.6) 



where = D^^ and 



^^={-1) -(2?+lK^|:+r--^(?' + 2?-47rV-8x+). (8.7) 
To show this, we apply il+ to (8.3). We have ^l+Ai(j)i^q{g) = •^i^+4>i,q{§>)- Thus, by (7.4), 

{1 1 ■ 1 ■ 

- 2^^^^ + ^'^rnre"^ cot 9 + -nre'-'^de - Trq- 



2sme 

+ Ir'd', - \mrdr + ^m^ - \rdr + ^m|e--^-^«'/'^L(rX,,(v[i^]). (8.8) 

On the other hand, invoking the first line of (7.12) and (7.14), we have 

(mcot 9 + do- ^) $L,,(v[i^]) = i{l-m + lX_i,,(v[i^]). (8.9) 

In the last two identities we set g = a[r], and note that ^(1) = 5m,q- Then we get the first 
identity of (8.6). In much the same way we get the second as well. 

Returning to the integral in (8.5), we see that it is equal to 

^°°<-i(r)<:;M-- (8.10) 



l-q + lJo ' ' r 

This procedure is valid only if Vq(r) tends to sufficiently fast as r tends to either or oo, 
which is in fact implied by the second line of (8.4). Hence 



dr 
r 



On the other hand we have, by the third line of (8.4) and by a formula of N.J. Sonine ([40, 
eq. (2) on p. 434]), 
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which gives, via either [14, eq. 4 of Section 6.576] or [30, (2.6.11)], 

and via (7.10), (8.5) and (8.11) we end the proof of (8.1). 

We turn to the surjectivity assertion. Thus, let F{u), tt e C^, be smooth and compactly- 
supported, and such that 

/ F{u)X(l)i,q{u)d''u = 0, for all (g,0 with IpI < |g| < L (8.14) 

We are to show F = 0. In view of (8.2) we may assume that F is radial, i.e., the Fourier 
expansion of F in u/\u\ has only one term, say, the g-th. With an obvious change in F, we 



consider instead 



/ , 

^0 



F{r)vUr)dr = 0, for aU {q, I) with \p\ < I, \q\ < I. (8.15) 



We then invoke that in [8, Lemma 5.1] more than (8.12) is proved; thus there are non-zero 



rf Jj; ly) such that 



l-inax{\p\,\q\} 

= Yl <,(j»r-"'^^^l^l'l^l>+'+^'^.+max{|p|,kl}-|P+.l+i(2^0- (8.16) 

Namely, we are given, for all integers I > max{|]3|, \q\}, 

F(r)r'+iK,+^_|p+,|(27rr)dr = 0. (8.17) 



/ 

^0 



We replace the Bessel factor by a well-known integral representation (see e.g., [30, (1.1.17)]), 
and find, after some rearrangement, that (8.17) is the same as 



^°°exp(-7rOr+'"'^+'''"'|^ F(r)r^-'^+|P+«lexp(-7rrV0^^}'^^ = 0- (§• 



18) 



Because of the completeness of polynomials over [0, oo), the member inside the braces should 
vanish for any ^ > 0; that is, for Re^ > by analytic continuation. Hence, the choice 
^ = 1/(1 + it), t e M, yields that the Fourier transform of a multiple of F{^) by a non-zero 
factor vanishes constantly. This ends the proof of Lemma 3. 

We now move to the proof of Lemma 4. It should be noted that with ordinary bounds 
for Bessel functions one may verify absolute convergence and analytic continuation needed to 
carry out the reasoning below. 

Thus, we consider the integral 

/•oo 

^iM= / v[{ry^'-^Ur. (8.19) 
Jo 
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The third hne of (8.4) gives 

Tl,,is) = ^l+^-2.(_l)min(0,p+,).-p-, + i(b + q\ - ^)) r^^^^ 



T{l-s + ^{\p + q\ + u)y 



with 



+ 



since for m e Z and — ^Iml < Res < ^ 



^0 



Jm{r)ir/2f'-'dr={-l) 



r(s + i|m|) 



(8.20) 



(8.21) 



r(l -s+ i|m|) 
We have the functional equation 

L,,,(s) = (-l)^-a,,_,(l-s), 



(8.22) 



(8.23) 



which is a result of the change of variable v v~'^ in (8.21). The necessary absolute conver- 
gence, and the meromorphic continuation to C of g(s) can be confirmed readily. Hence we 
have the local functional equation 



(-i)'-^ri,_,(s) = 7r2-^«(-i)™--(ipi'i'?i)r,,,(i - s) 

r(s + i(b + + iy))r{s + i{\p-q\-iy)) 



X 



r(l - s + ^{\p + q\+ u))T{l -s + ^{\p-q\- u)) 

(see [17, Theorem 6.4]). 

Then we observe, by convolving (8.22), that 



(8.24) 



POO 

/ A2^-V|p+,|(rA)J|p_,|(r/A)dA 

^0 



< — > 2 



s-3 



+ l{\p + q\ + u))T{ls + l{\p-q\-u)) 
r(l -ls + U\p + q\ + u))T{l -ls + U\p-q\- ^)) 



(8.25) 



is a Mellin pair, provided 2|Rez/| < Res < 1 — 2|Rez/|. Thus, denoting the left side by 
Ki/,p{r,q), we get, by (8.19), (8.24) and the Mellin-Parseval formula. 



(- 



/•oo 

■l)^-n-\LjA2) = 87r2(-l)"^^"(|p|'l^l) / K^^p{2T^\r,q)v[{r'^)rdr. (8.26) 

In this we set A = l^j with tt e C^, and multiply both sides by the factor (tt/lwD^^/H^p^^Hi^. 
On noting that $p,g(k[-/3, a]) = (-l)'-'^$p _q(k[Q;, /?]) or R^(j)i,q = (-l)'-«0i _q, we see by 
the definition (8.2) that (8.26) is identical to 

\u\-^XR^(l)i^q{u'') 

=47r(-l)"^^"(|p|'l«l) / K^,p{2TT\uvlq){uv/\uv\f'i%ct)i,q{v^)\v\'^d''v. 



=471 



<\p\^\^\)K^^p{27:\uv\,m) 



uv 



\uv\ 



2m 



%(t>i,q{v')\v\''d''v. (8.27) 
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We then invoke that Graf's addition theorem ([40, eq. (1) on p. 359]) gives, for any Z,z > 0, 



2mie 



E (-l)"^'^^'"'''"^'^<^|m+p|(^)J|^-p|(;^)e^ 

— oo 

= (-irJ.,(|iJe«' + .e-«|)(|ly±^) . (8.28) 



We apply this to the member inside the braces of (8.27), and find that the proof of (7.25) 
with (j) = (pi^q has been reduced to that of 

= (-ir^°°A-'J.,(2.M|Ae« + (Ae«)-'|) (^^±^^y dX, (8.29) 

with u — |w|e*^, which is, however, the same as [8, Theorem 12.1]. We end the proof; the 
extension to any smooth (p is immediate. 

It does not seem that the identity (8.29) had been tabulated before [8], a fact somewhat 
bizarre against its classical appearance. This integral representation of the Bessel function 
of representations of PSL2(C) is quite important, for it allows us to deal with test functions 
which not necessarily decay exponentially. This merit of (8.29) is indeed exploited fully in 
the proof of the spectral decomposition (7.28). The proof in [8] of (8.29) is conceptually 
involved, depending for instance on the Goodman- Wallach operator ([13]). The procedure 
above indicates the existence of a simpler approach, and in fact an alternative proof has been 
obtained in [7] [32]. Any extension of (8.29) is highly desirable. 

Lemmas 3 and 4 allow us to carry over the method of [11] to the complex situation, 
so that the proof of the Kloosterman-Spectral sum formula established in [8, Theorem 13.1] 
can now be proved in a simpler manner, although we have not worked out the details yet. 
Finally, we should mention that our argument seems to extend to Lie grpoups of real rank 
one; thus, the assertions due to R. Miatello and N.R. Wallach [22] are hoped to be included 
in our future discussion. 
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